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ABSTRACT 

We present a variational formalism for describing the dynamical evolution of an os- 
cillating star with a point-mass companion in the linear, non-relativistic regime. This 
includes both the excitation of normal modes and the back-reaction of the modes on 
the orbit. The general formalism for arbitrary fluid configurations is presented, and 
then specialized to a homentropic potential flow. Our formalism explicitly identifies 
and conserves both energy and angular momentum. We also consider corrections to 
the orbit up to 7/2 post-Newtonian order. 
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1 INTRODUCTION 

When a star orbits a companion, there can be a resonant ex- 
citation of stellar oscillations if the orbital period becomes 
close to a multiple of the oscillation period. This can hap- 
pen as a result of gravitational radiation or interaction with 
a third body, for example. This is a deceptively complex dy- 
namical system because there will be a back-reaction of the 
oscillation on the orbit. In this paper we consider a varia- 
tional approach to treating this problem. 

Variational principles for linear, adiabatic, non-radial 
stella r oscillations have been considered by several peopl e 
(e.g. IChandrasekharl |l964 iLynden-Bell & Ostrikeil h967|). 



Subsequently, a variational approach to the excitation of 



tides i n binar y systems was considered by Glingold & Mon- 
ag iianl (|l98C| ). However, their description was limited to 



polytropic equations of state, and did not account for or- 
bital evolution due to general relativistic effects which are 
significant for the long-term evolution of systems involving 
compact stars. In this paper, we consider a more general 
variational approach to tidal excitation. Although we limit 
our detailed discussion to irrotational flow, the formalism 
is presented in a way so as to allow generalisation to more 
complicated flows. 

In §2, we provide an overview of one variational ap- 
proach to the dynamics of perfect fluids. We then derive 
the most general Lagrangian for tidal excitation, and con- 
sider the special case of homentropic, irrotational flow. In 
§3, we derive the equations for tidal excitation in terms of 
the normal mode amplitudes, and obtain expressions for the 
conserved energy and angular momentum. We also consider 
corrections to the orbit due to general relativistic effects. In 
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§4, we discuss the applicability of our formalism. Finally, in 
§5, we reprise our conclusions. 



2 THE LAGRANGIAN 

2.1 Overview of variational fluid mechanics 

In this section, we outline some of the important results 
from a variational formulation of fluid mechanics. Our dis- 



cussion follows the review by Salmon (1988) where a more 



detailed exposition may be found. An important difference 
is that we have extended the formalism to accommodate a 
self-gravitating fluid in a non-inertial reference frame. 

The simplest variational formulation is to use a contin- 
uum version of the Lagrangian from classical mechanics. In 
this approach, the Lagrangian for the fluid is just the classi- 
cal Lagrangian for a system of particles distributed contin- 
uously in space. Let x(a, r) be the position, relative to the 
centre-of-mass, of the fluid particle identified by the label- 
ing coordinates a at time r. We shall distinguish between 
the time coordinates r and t. These are equal in magnitude, 
but partial derivatives with respect to r are at constant a, 
whereas those with respect to t are at constant x — in other 
words, d/dr corresponds to a convective derivative. There is 
considerable freedom in the choice of labeling coordinates, 
but it is convenient to choose them so that they are related 
to the mass density of the fluid by 



P = 



g(a) 
9(x)- 



(1) 



It should be noted that this just corresponds to the choice of 
a constant mass for the fluid particles (i.e. the mass density 
of the fluid is directly proportional to the number density 
of particles). This has the advantage that mass conservation 
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is implicit in our choice of a, as can be verified by a direct 
application of d/dr to (hi): 



dp 

dr 



+ pV ■ u = 0, 



(2) 



where u = 9x/9r, V is the gradient operator in x-space, 
and we have used the fact that the inverse of a matrix A 
may be written as 



ain|[^|[ 
dA„ 



(3) 



where ||^|| = det(A). We can now write down the La- 
grangian for the fluid as 



L, = da 



1 rm^ 9x 

2 I 9r ^ dr 



-i?('i4,S(a)) -$(x) 



a(a) 



(4) 



where R* is the location of the centre-of-mass, $ is the po- 
tential for external forces, and E is the specific internal en- 
ergy which is a prescribed function of the specific volume 
and the specific entropy S. Note that S depends only on the 
labeling coordinates a. This is, in essence, the perfect-fiuid 
approximation: 



dS_ 

dr 



0. 



It can be shown that the variation with respect to x of 
yields the Euler equation (see Appendix ^ . 

For cornputational purposes, it is more convenient 
to rewrite m) in Eulerian form. This is straightforward 
to accomplish by noting that the time-dependent map 
x — x(a, r) uniquely determines the inverse map a = a(x, t). 
Therefore, the requirement that the action be stationary un- 
der arbitrary variations (5x in the forward map is equivalent 
to the requirement that the action be stationary under vari- 
ations Sa in the inverse map. After dropping two total time 
derivatives, we can now write the fiuid Lagrangian as 



+ dx {p 



-u.u-iJ(p,5(a))-$(x) 



(5) 



where M* is the total mass of the fiuid. However, before 
we can consider the variation with respect to a of (^, we 
must express the velocity u as a function of a and its deriva- 
tives. Alternatively, we can include the relevant relations as 
constraints in the Lagrangian and then vary u and a inde- 
pendently. The required relations are given by 



„ da da , „. 



(6) 



which are the so-called Lin constraints. We may also include 
mass conservation d2|) as a constraint in the Eulerian form 



(H) of the Lagrangian. This gives us 



+ / dx i p 



iu.u-iJ(p,5(a))-$(x)-c |j 



(7) 



dp 
dt 



+ V ■ (pu) 



where ^ and are Lagrange multipliers, and we now consider 
the independent variations Su, Sa, S<^, Sp, and S(j>- 
The u variation of (^) yields 



u = CiVoi + V(l>, 



(8) 



which can be used to eliminate u from the Lagrangian. Af- 
ter dropping a time derivative and integrating one term by 
parts, (H) becomes 



dx 



^ da d6 1 
+ E{p,S{a)) + ^ 



(9) 



dp 
dt 



X R, 



where u is now just an abbreviation for (|8|), and the inde- 
pendent variations Sa, SC,, Sp and S4> '<^^^ to be considered. 

So far we have neglected the effects of self-gravitation, 
considering $ to be an externally imposed potential. We 
may now incorporate self-gravity in our formalism by in- 
cluding the Poisson equation. Our most general perfect-fiuid 
Lagrangian is then 



L. = -Af,R, 



dx 



^ 9a dS 1 ^ , r,, ^^ 

C- — + 7^+ oU-u + ^(P>5(a)) (10) 



dt 
* + $ 



dt 
1 

8^ 



:VVE' ■ V* 



dp 
dt 



X R, 



where ^ is the self-gravitational potential, and the indepen- 
dent variations are Sa, SC,, Sp, Sep, and S'^i. 

2.2 Coupling to a point-mass 

Consider now the problem in which the external potential 
$ is due to the gravitational field of a point- mass Mh- The 
Lagrangian for the whole system is therefore 

1 



-jd. 



. da dd> 1 ^, r^, ^-^ 

^ + ^ + - u + E {p,S{a)) + ^ (11) 



GMh 



|x-R| 

where R = R;, — R«. In the centre-of-mass frame of the 
system, (|ll[) becomes 



L = -pR^ 
- /dx 



C-^ + ^ + \^-^ + E{p,S{a)) + ^ (12) 



R 



dt 
1 
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where M = Mh + M, is the total mass, and fj, = MhM,/M 
is the reduced mass. 



2.3 Homentropic potential flow 

The variation of ( |l^ with respect to a gives us 

dc 



dE 
'dS 



dS_ 
da 



(13) 



We shall now derive the conditions on ^ for a homentropic 
potential flow. By definition, the velocity field for a poten- 
tial flow has the representation u = Vx' for some arbitrary 
scalar potential x' ■ We now define a new potential x such 
that x' ~ X + 4'- The velocity field can then be written as 

u = Vx + V<;!> = Va, + V0. 
oai 

Comparing this expression with we find that for a po- 
tential fiow 

. _ dx 
oai 

We now note that for a homentropic fluid it follows from ( p^ 
that x(a)''") = Xi{^) + X2(t). Therefore, after substituting 
(p^, we find that the Lagrangian ( p^ takes the form 



(14) 



dixi + 0) 



dt 



GMh 



R| 



8-kG 



-V* 



Since xi a^nd (fi only appear as the combination xi + ' 
can re-define xi + </> ^ <^ to obtain 













Ix - 


-Rl 



SttG 



-V* 



£;(p) + * 

V* H X ■ 

M dt 



(15) 



R 



It follows from our method of construction that all home- 
ntropic potential fiows can be derived from ( |l5| ) with as 
the velocity potential. 

We now note that the conservative nature of gravita- 
tional forces guarantees that, in the absence of viscosity, 
tidal excitation will not generate vorticity. Therefore, if the 
fluid starts out with zero vorticity then the flow will al- 
ways be irrotational. We can therefore use (^^ for problems 
involving non-rotating, homentropic stars (e.g. cold white 
dwarfs). From here onwards we shall assume this to be the 
case. 



3 EQUATIONS OF MOTION AND 
CONSERVATION LAWS 

3.1 The equations of motion 

We consider the first-order perturbations to (p, and p: 
0(x,t) = cj>o{r) + (l)i{x,t), 
*(x,t) = *o(r)-H*i(x,i), 
p{x,t) = po{r) +pi{x,t), 



where r = |x|, and we have assumed a static, spherically 
symmetric unperturbed fiuid configuration. Retaining terms 
up to quadratic order in the perturbations, we can separate 
the Lagrangian as 



L — Lcq -Lint + Lpc 



(16) 



where Leq describes the equilibrium (unperturbed) config- 
uration, Lpert describes the structure of the perturbations, 
and Lint describes the orbit-perturbation interaction. After 
making the perturbative expansion, we obtain 



R 



2' 
dx 

Lint = - / dx .J po (^01 + *i 

+ pi {ho + I'd 

+ 



GMh 
|x-R| 

1 d*o 9*1 Mh 



(17) 



(18) 



A-kG dr dr ^ M ^'"^ ' ^ 



+ ^PoV(?!)i ■ V(Pi + ^V*i ■ V*i 
2 ottG 



(19) 



where R = |R1, and ho and Cs are the unperturbed spe- 
cific enthalpy and adiabatic sound speed, respectively. The 
variations of Lcq with respect to R, po, and *o yield 



/lo + *o = 0, 

— -1- r — — = 4nGpo, 
dr \ dr I 



(20) 
(21) 

(22) 



where R is a unit vector. Together with an equation of state 
for the fiuid, (po|)"(p^ determine the unperturbed configu- 
ration. 

To determine the equations for the perturbations, we 
consider the variations of L with respect to 01, pi and "ti. 
These give us 

Pi-l-V-(poV(;!>i) = 0, 

„2 



(pl + —pi + *1 = ^ + —X • R, 

po - R| M 

V^*i = AttGpi. 



(23) 
(24) 
(25) 



It is straightforward to show that, with M^ — 0, (|23|)-(|2^) 
are the conventional equations for the normal modes of a 
non-rotating, homentropic star (see Appendix^. 

The equation for the orbit, including the back-reaction 
of the perturbations, is obtained by considering the variation 
of L with respect to R: 

where we have used the fact that 

dx /9ix = 

(this is equivalent to choosing the origin of the coordinates 
X to be the centre-of-mass of the fluid). 
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3.2 Amplitude formulation 



3.3 Energy and angular momentum 



It is convenient to write the equations for tidal excitation in 
terms of the normal mode amplitudes. Let ^ be the physi- 
cal displacement field of fiuid elements within the star. Inte- 
grating ( p^ ) with respect to time and setting the integration 
constant to zero, we get 



Pi = -V- (po4). 

Taking the gradient of (p^, we have 

i' + v 



—pi + *i = GMhV 
Po / V |x 



R| 



(27) 



(28) 



The second term on the left side of 
linear, spatial operator on ^: 

2 



28h can be written as a 



^^V ■ (poO 
Po 



G / dx 



(29) 



where we have used (^^ and ( psj ) to write pi and ^i in 
terms of ^. We can solve for the eigenfunctions ^j(x): 

where are the eigenvalues. It can be pro ved that the op- 
erator X> is Her mitian with respect to mass ( Chandrasekhar 



1964; Cox 198C). We are therefore guaranteed that its eigen- 



values are real, and it is generally assumed that its eigen- 
functions form a complete set and can be made mutually 
orthogonal. Hence, with the appropriate choice of normal- 
ization, 



Expanding ^ as 



dx poCf ■ = Sjf 



|(x,t) = ^^(t)|,(x) 



(30) 



(31) 



in (ESI) and then projecting out a single mode gives 



A, + Lo^Aj = /, (R) -f . / dx po^; , (32) 



where 



/j (R) = GMh /dxpoS; ■ V 



|x-R|' ^''^ 

(see Appendix ^). This is just a forced, harmonic oscilla- 
tor with frequency Uj. The second term on the right side 
of ( ^ ) is non-zero only for for monopolar (Z = 0) modes, 
and exactly cancels the first term for that case. This is a 
mathematical statement of the fact that monopolar modes 
are not tidally excited. 

The equation for the orbit (^), after substituting ( |27| ) 
and then performing an integration by parts, becomes 



Expanding { in terms of the normal modes, we find 

■A 
'aR 



To find the conserved energy £, we calculate the time-time 
component of the energy-momentum tensor from (^ij) using 



dC 



d{diqk 



(35) 



where £ is the Lagrangian density, and are the generalized 
fields. A straightforward evaluation gives 



dxT* 



= 2^^ 



R 



GMh / dx 



Pi 



ix-R| 



(36) 



+ \ [ dyi ( poV0i ■ V01 + —pi + pi*i 
2 J \ Po 



where we have used ( p5| ) to eliminate the gravitational self- 
energy of the perturbations. The total energy is the sum of 
three components: orbital, perturbation, and coupling. The 
various pieces are easily identified in (p6[). 

We can, without loss of generality, assume the orbit to 
be in the equatorial plane. The conserved angular momen- 
tum Lz can then be calculated from the appropriate com- 
ponent of the energy-momentum tensor as 



= / dx 



, a<pi 

dx pi^, 



(37) 



where (p and u are the azimuthal coordinates associated with 
X and R, respectively. 

It may be noted that the canonical form (^^ of the 
energy-momentum tensor is not manifestly symmetric. How- 
ever, it is well-known that the tensor can be made symmetric 
by the addition of a suitable divergence term: 



ikj 



We do not need to do this since integral quantities such 
as (Is^) and (^tI) a re unaffected by such a transformation 
(Landau & Lifshitz 1975). 

The energy and angular momentum associated with the 
normal modes also take on relatively simple forms in terms 
of the time-dependent amplitudes Aj. From (^6|), we know 
that the energy associated with perturbations is 



^pcrt — 



dx 



pol ■ ^ + Pi —pi + *1 



po 



(38) 



Substituting (^^ into (^^ and then integrating the second 
term in the integrand by parts gives 



fpcrt = - / dx Po 



|.| + |.V ^pi-h*i 
po 



Note that the second term in the integrand now involves 
the same linear operator that we used to define the normal 
modes. Expanding ^ as in (^l|) and using the orthonormality 
relation (pOl), we find that the energy associated with mode 
j is just 



A, 



(39) 
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From (|3^) and (|27[), we know that the angular momentum 
associated with perturbations is 



dip 



Performing an integration by parts, we get 



- pert 



d(fi 



(40) 



Once again expanding ^ as in (|31 
relation (pOl), and the fact that 



using the orthonormality 



dip 



(see Appendix^, we find that the angular momentum as- 
sociated with mode j is just 



-imAjAj . 



(41) 



We shall now derive a simple relation between the energy 
and angular momentum associated with an isolated mode 
(i.e. with fj =0). In that case, from (^) we have Aj{t) oc 
e*"j*. Therefore, from ( psj ) and (|l|) we get 



mujj\Aj 



(42) 



This relation is to be expected on physical grounds as fol- 
lows. If we consider tidally exciting the mode at resonance in 
a circular orbit, then the rate at which energy is transferred 
to the mode is just 



dt dt m' 

where r is the torque exerted by the perturbing mass, and 
Q. is the orbital frequency. Integrating this equation with 
respect to time and setting the initial mode energy and an- 
gular momentum to zero, we obtain ([43). 

Finally, the conserved energy (pGHand angular momen- 
tum (^), written in terms of the amplitudes, are 



£ 



-pie 



R 



and 



^imAjA]. 



(43) 



(44) 



3.4 Post-Newtonian corrections 

For problems involving compact binaries, general relativistic 
effects may be important. The corrections to the orbit are 
partly periodic, but there are important secular effects, such 
as periastron advance and gravitational radiation reaction, 
which must be taken into account for even a qualitatively ac- 
curate treatment of the problem. The most straightforward 
scheme, and the one we adopt, is to treat the corrections as 
perturbations to the Newtonian orbit, neglecting all finite- 
size effects. In this scheme, we simply add the appropriate 
corrections to the orbital acceleration, energy and angular 
momentum at each time-step. Detailed derivations of thes e 
corrections may be found elsewhere (cf. Iyer & Will 1995). 



We only list the relevant equations here for reference as they 
will be used in a forthcoming p ublication. The equations are 
taken from |lyer fc Will) ( [l995t ). 

The first- and second-order corrections to the orbital 
accelerations are 



R 



(1) 



GM f 

i?2 [ 



2(2 + 7,) 



GAI 



+ (1 + 37?) 



3 
2 



R-2(2-r?) 



(45) 



( 5 



d(2) 



R 



GM / 

' i?2 \ 



+ 7,(3-4,7)^ + ^7,(1-3,?)^ 



-r,(3 - 4r?) 



(2 + 25r/ + 2 J," 



,4 



1 ^ ,GMR2 

2^(13 -47;)^3^ 



GMR^ 



d^R 



(46) 



R 



R^ 



7,(15 + 47,)^ - (4 + 417, + 877") 



2^ GM 



where r] = p/M. These corrections are formally conservative 
in the sense that it is still possible to define a conserved en- 
ergy and a conserved angular momentum to second-order in 
ie /c^. The appropriate corrections to the Newtonian quan- 
tities are 



C-p^ — pL, 



3/, o ^R'' 1/n ^GMR2 
0(1-37,)— + -(3 + 7,)- 



c-^R 



1 GMR^ 1 fGM\' 
^ 2^^ c^R 2 I "^j 



(47) 



C-pjq pi- 



^(1-^7, + 137,^)^ 



+ -(21 - 237,- 277," 



GMR^ 



n Q . OMR" 



(48) 



+ -(14 - 557, + 47,^) ' ' 



1^ c3R J 



(49) 
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and 



r(i) 



r{2) 



1, ,11^ , ,GM 
-(1-3,)- + (3 + ,)^ 

2(7-10.-9,^)^3^ 



(50) 



(51) 



1 , ,GM7?2 
-2^7(2 + 5,)^ 

+ i(14-41. + 4,^)(||)^ 

+ |(l-7r, + 13r,^)^ 

At 5/2 post-Newtonian order, the first dissipative terms 
arise. Therefore, after second-order it is no longer possible 
to define a conserved energy or angular momentum. The 
5/2- and 7/2-order corrections to the orbital acceleration 
are given by 



R.RR 

where 



5' \ R'^ J \c^Rj 



R 

(A5/2+A7/2)-R. 

C 



+ (-B5/2 + -B7/2) — 



S5/2 = 6 



2GM 
GM 



R^ 



c^R 



15- 



.R^ 



87 



R* /5379 ^ 1^ ^ GMH 



V 28 ^ 3 ^ 



+ 



V 4 



-H3377 



35 
2 



,7?* /160 



55 



c*R 
GMFf 
c*R 
2 



(52) 

(53) 
(54) 

(55) 



B. 



7/2 



27R^ /4861 58 \ GMR^ 
3/,o ,7?^R^ /2591 „^ \ GMR" 



(56) 



.T?-* 1 /776 _ \ /GM\ 
^ + - h 55r) — — 



f(l-7,)^ + -^^+55,^^^,^^ 



4 APPLICABILITY 

Our formahsm, as developed in the previous section, is valid 
only for non-rotating, homentropic stars. However, the La- 
grangian (|l^) has more general applicability to any problem 
where dissipative eflects may be neglected. It is therefore 
possible, in principle, to generalise the formalism to include 
rotating fluid configurations and non-homentropic equations 
of state. In this section we shall discuss further limitations 
of the formalism developed in § ^ 

The two most important cases where our formalism is 
not applicable are regimes of strong non-linear eflects and 
gravity. In particular, one fundamental assumption upon 
which our mode analysis is based is that the underlying 
fluid configuration remains spherically symmetric. This is 
no longer the case, for example, when the star fills its Roche 



lobe and is tidally disrupted. We can formulate this as a 
quantitative limit using an approx i matio n to the eflective 
Roche lobe radius due to Eggleton (198J): 



R^ 
R 



< 



0.49 



2/3 



0.6g2/3 +in(l + gl/3)' 



(57) 



where q = M,/Mh- 

The limitation in the regime of strong gravitational 
fiel ds is due to the post-Newtonian corrections summarised 
in § 3.4 , In particular, the perturbation series in the order pa- 
rameter ti? /c? , represented by the corrections, fails to pro- 
vide a satisfactory approximation in the strong-field limit. 
The point at which the p ost-Newtonian cor rections fail is 
roughly at J? = 10 GM /c^ (iBrady ct ai]|l998|). Higher orde r 



corrections have been calculated ( Gopakumar et al. 1997 ), 
however they are not expected to significantly increase the 
region of validity of our formalism. To evolve the orbit in 
the strong-field regime, a fully relativistic treatment of the 
two-body problem is required. 



5 CONCLUSION 

In this paper, we have described and developed a variational 
approach for handling a specific dynamical problem — the 
oscillation-orbit interaction of a non-rotating star with a 
compact object companion. Our treatment is presented in 
some generality because we believe that this method can 
also be employed in other problems associated with accre- 
tion disks and extra-solar planets. We recover the standard, 
Newtonian equations of fiuid dynamics for irrotational fiow 
as wells as normal mode equations for small oscillations. The 
power of the Lagrangian approach is made manifest in the 
identification of algebraic expressions for the modal energy 
and angular momentum, and for the equations of motion. It 
also suggests an approach to handling more complex prob- 
lems, including those involving strong-field gravity and ro- 
tational fiow. 

In a forthcoming paper, we will use this approach to 
compute the dynamical evolution of a white dwarf in orbit 
about a compact object as it spirals inward under the action 
of gravitational radiation. 
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APPENDIX B: THE NORMAL MODES 

We expand pi, and |x — R|^^ in terms of spherical 
harmonics: 



l,m 
l,m 

*i(x,i) =X*,™(r,t)y,,„(x), 



Ix-R| 



APPENDIX A: THE EULER EQUATION 

We shall assume, for simplicity, that the centre-of-mass 
frame of the fluid is inertial. Then, in the centre-of-mass 
frame, the Lagrangian (W) becomes 



L, = da 



9x 

2 V 9r 



i5(Q,S(a))-$(x)], 



where a = p ^ is the specific volume. The variation with 
respect to x yields 



dE d 
da duj 
da 



da 



d{dxi/daj) 
d fdE 



dxi ' 



(Al) 



d{dxi/daj) daj \ da 
We can show that the second term on the left hand side of 



(Al) is zero as follows. From the definition of a, we know 
_ 9(x) 1 dxi dxj dxk 



that 



9(a) 



dm dam dan ' 



Arfken & Weber 



where egk is the three -dimensional Levi-Civita symbol (cf. 

1995). Taking the derivative of both sides 

dxj dxk 



with respect to dxi/daj, we get 



da 



d{dxi I daj 



' dam dan 



And, finally, taking a derivative with respect to aj gives us 



d 
daj 



da 



d{dxi/daj 



d Xj dxk _ Q 



daidam dan 



Next, using the deflnition of a and the identity (^), we note 
that 



da 



d 



d{dxi/daj) daj 



daj d 
dxi daj 



d 
dxi 



Therefore, (Al) now becomes 



d'x. 



dP 

dxi 



9$ 

dxi '' 



where P = —dE/da is the pressure. Recalling that d/dr 
corresponds to a convective derivative in the Eulerian de- 
scription, we see that this is just the Euler equation. 



Jackson 



(cf. 

grating over angles, (|23|)-(|2^) become 

1(1 + 1 



1999). After inserting the expansions and inte- 

PO(l>lm = 0, (Bl) 



19/2 d<f>lm 
Plm + K PO—^ 

dr \ dr 



km + —Plm + ^Im 
PO 



GMh 47r 
R 21 + 1 \R 



Ylmi'R-) 



+ '5n-^r|R|y4(R) 



J.2 Qj. 



dr 



1(1 + I) 



(B2) 



(B3) 



The flrst term on the right hand side of (B2) is a forcing 
term that couples the modes to the gravitational potential 
of the point-mass. The second term, which is only present 
for dipolar (/ = 1) modes, exactly cancels the first term for 
that case. We can therefore conclude that dipolar modes are 
not tidally excited. This was to be expected since the origin 
of the coordinates x is the centre-of-mass of the fluid. 



The temporal Fourier transforms of (B1)-(B3), with 
Mh = 0, give 



dr ^ dr 



2 d 1(1 + 1) ' 

-;:2— 



iuj(f>lm H -plm + ^im ~ 0, 

PO 



1 d 2 d 1(1 + 1) 
dr dr r'^ 



Plm = —IL^Pln 



4nGpim, 



(B4) 
(B5) 
(B6) 



where pim, 4'im and ^im are the temporal Fourier trans- 
forms of pim, 4'im, and 'i'lm, respectively. We can use (B5) 
to eliminate jiim from (B4) and (BC). This yields the two 
second-order equations 



d 2 d 1(1 + 1) LO^' 

-J Po-r — + — 



r^po dr dr 



l^'^lm, 



(B7) 



1 d 2 d 1(1 + 1) 4ttGpo 
r^ dr dr c? 



With the deflnitions 

1 d<j)ln 



-I 5 (plm- 

(B8) 



iur dr 
1 



- - ^ 1 

r}2 = -. <Plm, 

igr 



-4- 



gr 



1 d*!,, 



g dr 



© 2002 RAS, MNRAS 000, 1-| 



8 Yasser Rathore, Avery E. Broderick and Roger Blandford 



and 



U E 



d\nM 
d\nr ' 
dlnPo 
d In po ' 

-UJ , 



V = - 



dlnPr 







dlnr 



where 



M{r) 

air) = 



= / dr' 4:Tvr' po{r'), 

Jo 
GM{r) 



and after some manipulation, (B7) and (B8) give the four 
first-order equations 



d^ 

dr 
^ dfj2 

dr 
^ dfjs 

dr 
^ drjA 

dr 



V_ 



■ (1 - U)f,-A + 774, 

1(1 + 1) - 



uv _ 
-T^m + 
i 1 



uv 



V2 + ^V3, (B9) 
i 1 



(BIO) 
(Bll) 
(B12) 



These are the conventional equations for the normal mode 
structure of a non-rotatinR, h omentropic star (cf. Cox 
Kippenhahn fc Weigert||l99C|). 



1980 



To have a well-posed problem, we obviously need to 
specify four boundary conditions. At the center of the star, 
we require that the variables 77; be well-behaved. Expanding 
in a power series around r = 0, we have 



(|B13|)-(|B16[), (|B9|)-(|B12|) consitute an eigenvalue problem 
for the normal modes. 

We now turn to the problem of determining the phys- 
ical displacement of fluid elements from the unperturbed 
configuration in terms of the r^i. Expanding (the tempo- 
ral Fourier transform of the perturbation to (p) in terms of 
the normal modes, we have 

= ^ Anlm{(^)(i>uhn{r)Yl,„{x.), 

where 4'nim is a normalized eigenfunction, and we use the 
subscript n to distinguish between the various modes cor- 
responding to the same l,m. Using the fact that ^ = V<?!>i, 
where ^ is the physical displacement of fluid elements, we 
find 

('^)in!m(x), (B17) 

n,l.m 

where 

|„i„(x) = ^rfiix + ^^772vj yim(x). 

The displacement fiel d $( x,t) is just the temporal inverse 
Fourier transform of (B17). 

, it is straight- 



Note that by taking the gradient of (| 

forward to see that the normal modes ^nim satisfy the op- 
erator equation 



V ( — pi + *i 

po 



) a 

'r . 



Substituting into (B£)-(pi2|), and using the fact that 



lim [7 = 3, 

1 — .0 

lim V = Q, 

1 — ,0 

lim C = constant, 

1 — »o 

we find that the only non-vanishing coefficients correspond 
to a = / — 2, and 



(B13) 
(B14) 



which constitute our boundary conditions at the center. At 
the surface, we require that ^im satisfy the Laplace equa- 
tion. This gives 



(B15) 



at r = r*, as our third boundary condition. Finally, from 
( |B9[ ) and using the condition that 

GAf, 

Po 



lim V 

r — 



hm = OQ 



(cf. Cox 198C ) , we get a fourth boundary condition that 

V2 = fii + m, (B16) 
at the surface. Together with the boundary conditions 



APPENDIX C: THE OVERLAP INTEGRAL 

Integrating ( [33[ ) by parts and substituting ([27[), we find 



From 



/, (R) =GMh / dx 



we know that 



|x-R| 



*;(R)^-G/dx^ 



Therefore, we have 



/,(R) = -Mft**(R). 



(CI) 



It is straightforward to write 'I'j(R) in terms of the dimen- 
sionless variable 7)3 defined in Appendix ^ using the solution 
to the Laplace equation in spherical coordinates: 

*.(R) = ^'?3(p,) 

Therefore, the overlap integral ( |Cl[ ) is 

/.(R) - ~^^UR*) (f )'^Vr„(R). 

The above expression for the overlap integral explicitly 
demonstrates that in the absence of a perturbation to the 
gravitational field there is no mode excitation. This is, of 
course, a direct consequence of Newton's third law. Assum- 
ing the orbit to be in the equatorial plane, the expression 
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becomes 

GMhM, _ 

i+i (C2) 



X I — I e™", 



where u is the azimuthal coordinate associated with R. 

This paper has been typeset from a Tj^X/ I^TjEX file prepared 
by the author. 
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